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Abstract. Using a large sample of Mgll absorbers with 0.4 < z < 2.2 detected by Nestor et 
al. (2005) in the Early Data Release of the SDSS, we present new constraints on the physical 
properties of these systems based on two statistical analyses: (i) By computing the ratio between 
the composite spectra of quasars with and without absorbers, we measure the reddening effects 
induced by these intervening systems; and (ii) by stacking SDSS images centered on quasars 
with strong Mgll absorption lines and isolating the excess light around the PSF, we measure the 
mean luminosity and colors of the absorbing galaxies. This statistical approach does not require 
any spectroscopic follow up and allows us to constrain the photometric properties of absorber 
systems. 
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1. Introduction 

Since it was realized that strong metal absorption lines detected in quasar spectra are 
due to intervening galaxies (e.g. Bahcall & Spitzer 1969, Boksenberg & Sargent 1978, 
Bergeron & Boisse 1991), significant progress has been made in our understanding of 
absorber systems. A number of studies have constrained some of their physical properties 
such as temperature, abundances and kinematics (for a review see for example Churchill’s 
contribution). However, the precise nature of absorber systems is still elusive and more 
observational constraints are needed in order to improve our understanding on how these 
systems trace various structures such as galactic disks, halos, winds, etc. This would allow 
us to include and use absorbers in our global picture of galaxy formation and evolution. 

In this contribution we illustrate how imaging and spectroscopic data from the Sloan 
Digital Sky Survey (SDSS) can be used to constrain the nature of absorber systems. 
Having selected a population of quasars with strong absorbers as well as a population 
of reference quasars (§2), we construct their composite spectra and isolate the redden¬ 
ing effects due to the presence of absorbers (§3). In §4, we show how the photometric 
properties of absorbing galaxies can be constrained by measuring the statistical excess 
of light present around quasars in SDSS images. We summarize the results in §5. 

2. The data 

2.1. The Mgll absorber sample 

We use the sample of Mgll absorption line systems compiled by Nestor et al. (2005) 
based on SDSS EDR data (Stoughton et al. 2002). Here, we focus our analysis on strong 
systems with a rest equivalent width Wo(A2796) > 0.8 A. In this range, the detec¬ 
tion completeness is very high and multiple systems along the same line-of-sight are 
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Figure 1 . Left: redshift distributions of the Mgll absorbers with W(A2796)o > 0.8 A (solid line) 
and background quasars (dashed line). Right: rest equivalent width distribution of the selected 
absorbers. 


rare. In addition, we focus our reddening measurements on the strongest systems, with 
W(A2796)o > 2.5 A. In this section we briefly summarise the main steps involved in the 
absorption line detection procedure. For details on the quasar and absorber catalogs, we 
refer the reader to Schneider et al. (2002) and Nestor et al. (2005). 

The SDSS EDR provides 3814 QSO spectra, approximately 3700 of which are of QSOs 
with sufficiently high redshift to allow the detection of intervening Mgll absorption lines. 
The continuum-normalized SDSS QSO spectra were searched for Mgll AA2796,2803 dou¬ 
blets using an optimal extraction method employing a Gaussian line-profile to measure 
each rest equivalent width Wq. The identification of Mgll doublets required the detection 
of the 2796 line at more than 5cr and one additional line detected at least at 3cr. Only 
systems 0.1c blueward of the quasar redshift and red ward of Lya emission were accepted. 
Systems with separations of less than 500 km/s were considered as a single system. The 
final sample comprises 813 QSOs with Mgll absorbers. We show the redshift distribution 
of the absorbers with W(A2796)o > 0.8 A (solid line) and that of the background quasars 
(dashed line) in the left panel of Figure|T| The right panel shows the rest equivalent width 
distribution Wo of the 2796 A line. Among these systems, 55 have W(A2796)o > 2.5 A. 

2.2. The reference quasar sample 

In order to isolate the effects induced by the presence of absorber systems we use a 
control sample made of quasars without absorbers. To avoid selection biases, the two 
quasar populations must have the same redshift distribution and the same absorption¬ 
line detectability. The main steps of our object selection are as follows: for each quasar 
with an absorber system, we randomly look for a quasar without any detected absorber 
(down to the limiting value of the line finder, i.e. W(A2796)o = 0.3 A ) such that their 
redsliifts are similar within 0.1. In addition, we require that this quasar has a high 
enough S/N to allow the detection of absorption lines of similar strength as the absorber 
system at wavelengths corresponding to z a bs ■ Thus our procedure ensures that the lack of 
absorption at z a bs in the reference QSO spectrum is real. The quasars selected in this way 
are called reference quasars in the following. Contamination by undetected absorbers at 
^ A Zabs is expected to statistically affect both samples in the same way. We match each 
quasar with absorbers to two reference quasars, which allows us to decrease the noise 
level. 
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Figure 2. Composite spectra of quasars with and without strong Mgll absorbers with 
Wo > 2.5A, in the observer (left) and absorber (right) rest-frame. The lower panels show the 
ratio of the composite spectra. The fitted curve used in the lower panels corresponds to the 
SMC reddening curve. 


3. Reddening properties from composite spectra 

In this section we investigate the reddening and extinction properties induced by our 
selected population of absorbers to their background quasars. For each sample of quasars 
defined above we create composite spectra in the observer- and absorber rest-frames. 
The composite spectra are obtained by computing the geometric mean of all spectra 
renormalized so that the global signal is not dominated by the brightest objects and 
taking into account the masks provided by the SDSS pipeline. For computing the com¬ 
posite spectrum in the absorber rest frame, sky lines are masked and the wavelength 
shifts applied to the quasars with absorbers are also used for the reference quasars. The 
results are shown in the upper panels of figure 0 The composite spectra are shown in 
the observer and absorber rest-frame in the left and right panels. The corresponding 
spectra have been normalised to the same flux at A = 8500 and 3500 A, respectively. 
As can be seen in both panels, quasars with absorbers appear to be redder than their 
reference quasars. In order to quantify this difference, we have computed the spectrum 
ratio and displayed the result in the lower panels of the figure. In the absorber rest frame, 
this ratio provides a direct measurement of the mean extinction curve of our sample of 
absorber systems. In the observer rest-frame, this quantity has been convolved with the 
redsliift distribution of the absorber systems. As we can see, the 2175 A excess seen in 
the reddening curve of the Milky way or the LMC is not present. We find that the SMC 
reddening curve provides a good fit to the mean extinction curve of Mgll absorbers’). 
This is shown with the smooth curve. In each case, we have normalized the ratio such 
that it is unity at sufficiently large wavelengths. This allows us to estimate the expected 
amount of extinction as a function of wavelength. 

Using the SMC reddening curve as a fitting template, the measured reddening allows us 
to constrain the mean product (( k/ksAic) Np) where Usmc is the SMC dust-to-gas ratio 
and Nq is the hydrogen column density. From the results obtained in the observer and 


t We use the template provided by B. Draine at http://www.astro.princeton.edu/~draine 
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absorber rest-frames, we find (( k/ksMC ) Nh) = (2.59±0.32) x 10 21 and (2.08±0.52) x 10 21 
atoms cnr 2 , respectively. Note that the noise contributions are different in each case. 
Our result indicates that Mgll absorbers with IF(A2796)o > 2.5 A generally correspond 
to Damped Lya systems. We find the reddening values to strongly depend on the Mgll 
rest equivalent width. Detailed results will be presented in Menard et al. (2005). 

4. Photometric properties from image stacking 

In this section we show that producing a high signal-to-noise (S/N) average image 
of the quasars with absorbers allows us to measure the photometric properties of the 
absorbing galaxies. For the two samples of quasars defined in section 2, we stack the 
corresponding SDSS images centered on the quasars. This procedure consists of four 
basic steps: shifting to make images superposable, rescaling the intensity to a uniform 
photometric calibration, masking of all unwanted sources, and averaging. 

We use the raw “corrected frames” publicly available in the Third Data Release of the 
SDSS (Abazajian et al. 2005). The image shift is done using a bi-cubic spline interpolant 
and the photometric calibrations provided by the DR3 catalogs are used to rescale the 
intensities. We then generate masks blanking sources that are likely to be foregrounds 
with respect to the absorber. For each QSO we create three independent masks in g, r and 
z-band, that are eventually combined in a final mask. In each band, we first subtract from 
the corrected frame the Atlas Image of the QSO produced by the PHOTO pipeline (Lupton 
et al. 2001). As the Atlas Image contains all and only the pixels that are attributed to 
the QSO by the de-blender, we are left with a clean image where all sources but the 
QSO are present. On this image we run SExtractor (Bertin & Arnouts 1996) to find 
all candidate sources. We mask all objects with a flux exceeding that expected for an 
extreme case of a luminous blue galaxy: a stellar population with Mj < —22.0 produced 
by a burst 10 7 Gyr long, observed immediately after the burst end at the redshift of 
the absorber as predicted by Bruzual & Chariot’s population synthesis models (Bruzual 
& Chariot 2003). Our final results do not critically depend on this choice of different 
spectral energy distributions. Finally, the masks in the three bands are combined and 
conservative masks are applied to all objects classified as stars by PHOTO with r < 20.0. 
Non-saturated stars are masked out to their 25 mag arcsec -2 isophotal radius (r-band), 
while saturated stars are masked out to three times this radius. 

After subtracting the background value evaluated on the masked image, the shifted 
and intensity-calibrated images are averaged, excluding all masked pixels. Azimuthally 
averaged surface brightness profiles are extracted in a series of circular apertures which 
are optimized to integrate on larger areas at larger radii. A more accurate background 
level for these profile is recomputed on the stacked image. The complete covariance matrix 
of the extracted background-subtracted surface brightness profile is evaluated using the 
jackknife technique. 

In order to isolate the excess light due to intervening absorbers, we subtract the mean 
point spread function (PSF). The PSF depends on several parameters: the time of the 
observation, the position on a given camera column, the color of the object, etc. In order 
to build an effective PSF which is representative of the flux-averaged PSF of our QSOs 
extended to large angular distances, we stack stars selected to match each individual 
QSO. We first require stars to be observed in the same observing run and same camera 
column as the QSO. These stars are then ranked to minimize a combination of the 
following quantities: (i) the difference of Gaussian FWHM with respect to the QSO; (ii) 
the angular separation from the QSO; and (iii) the differences in colors with respect 
to those of the QSO. To map the PSF with reasonable S/N at large distances, bright 
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Figure 3. The statistical light present around quasars with (upper panels) and without Mgll 
absorbers (lower panels) measured in the g , r and i-bands after having subtracted the PSF. The 
light excess seen around quasars with absorbers corresponds to the light originating from the 
absorbing galaxies. 


unsaturated stars are also preferred (r < 17.0 mag is however required). For each QSO 
in each band, the star that best fulfills all these criteria is chosen. To obtain the correct 
flux-weighted average PSF, the intensity of each star is rescaled to the same intensity of 
the corresponding QSO before stacking. 

In FigO we show the surface brightness residuals measured in the g , r and i-band 
after the PSF subtraction for the sample of 810 QSOs with absorbers (upper panels) and 
for the reference QSO sample (lower panels). Note that three quasar images could not 
be used because of their proximity to a bright star. The stacked PSFs are normalized 
with respect to the QSO profiles within an aperture of 4 pixels (1.6"). The vertical scale 
is linear in flux intensity and the horizontal solid line is the background level. The left 
vertical axis displays the surface brightness in mag arcsec -2 (in excess or deficit). Error 
bars are computed from the jackknife covariance matrix and include the budget from the 
background uncertainty. Central normalization uncertainties are found to be negligible. 
The upper panels show a clear surface brightness excess. It is detected at 5.5, 7.1 and 
9.2cr in the g , r and i-band, respectively. Such an excess is not seen for the reference 
sample. In this case, the results indicate a marginal surface brightness excess within 4", 
possibly due to the flux contribution of the quasar host galaxy. Larger datasets will allow 
us to measure this effect more precisely. The excess light that cross-correlates with the 
presence of Mgll systems provides us with an estimation of the mean light contribution 
of the absorbing galaxies. In addition it should be noticed that the width of this light 
excess corresponds to the characteristic size of Mgll gas halos around galaxies. 

We have investigated how the signal depends on redshift. Within 10" the integrated 
i-band absorber magnitude range from i = 21.50t 0 ;gg at 0.37 < z ^ 0.82 to roughly two 
magnitudes fainter at 0.82 ^ z ^ 2.2, which is consistent with the expected dimming as 
a function of distance. More detailed results are presented in Zibetti et al. (2005). 

Focusing on the redshift range 0.37 < z < 0.82 we have investigated the photometric 
properties of the absorbing galaxies. Within a 10"aperture (where the S/N is higher) 
we find the average integrated magnitudes to be g = 23.48+q2o, r = 22.27+0.12 an d 
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i = 21.50lo 08 > an d the corresponding colors read g—r = 1.22±0.26 and r—i = 0.77±0.15, 
indicating a population of blue galaxies. More detailed results on the implied stellar 
population will be presented in a forthcoming paper. 


5. Conclusion 

Using a sample of over 800 Mgll systems with 0.4 < 2 < 2.2 we have presented two 
statistical analyses constraining the properties of these absorber systems. After having 
carefuly selected a population of reference quasars, we have measured the mean red¬ 
dening curve induced by these intervening systems by computing the ratio between the 
composute spectra of quasars with and without absorbers. We do not find any presence 
of the 2175 A bump present in the Milky way and LMC reddening curves and we find 
the SMC reddening curve to provide a good fit to the measured reddening. Assuming the 
mean Mgll absorbers to have the same metallicity and dust-to-gas ratio as the SMC, we 
measured the mean hydrogen column density of these systems to be ~ 10 21 atoms cm -2 . 

In addition, we have constrained the mean photometric properties of Mgll absorbing 
galaxies by measuring the statistical light excess present around quasars with absorbers. 
This signal is detected at 5.5, 7.1 and 9.2 a in the g , r and i-band, respectively and such 
an excess is not seen for the reference sample. This light contribution behaves as expected 
as a function of scale and redshift, and the measured colors indicate that the mean Mgll 
absorbing galaxies is a late-type galaxy. 

The results presented in this contribution are based on EDR data. Future analyses 
will provide much more accurate constraints and will allow us to investigate the absorber 
properties as a function of redshift and for different ion species. 
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